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Introduction 


The  major  accomplishment  reported  herein  is  the  invention  of  a 
two-beam  accelerator,  named  "Twobetron"  [1]  that  was  supported  by  this 
contract.  A  patent  is  being  granted  [2].  This  accelerator  uses  the  power 
of  an  intense  annular  relativistic  electron  beam  to  drive  a  low  current  on- 
axis  pencil  beam  to  high  energies  for  industrial  and  medical  applications. 
For  the  purposes  of  this  report,  high  energy  is  in  the  1  to  1 0  MeV  range 
rather  than  the  GeV  range  in  the  high  energy  physics  community. 

Summarized  below  are  the  concepts,  theories,  simulation  results, 
studies  of  beam  modulation,  and  initial  experimental  studies.  The  details 
are  given  in  the  references.  The  ones  supported  by  this  contract  are 
attached  to  the  end  of  this  Report. 

Theory 

Concept 

In  the  Twobetron,  an  intense  annular  relativistic  electron  beam  with 
current  Id  travels  down  a  series  of  pillbox  cavities  at  radius  r©  to 
accelerate  a  pencil  beam  with  current  Is  «ld-  Both  the  driver  beam  (Id) 
and  the  secondary  beam  (Is)  3re  modulated  at  a  frequency  wq,  the 
frequency  of  the  TMo20  mode  of  the  cavity  [Fig.  1].  The  driver  beam  is 

always  decelerated  by  the  mode,  by  conservation  of  energy.  Since  the 
electric  field  of  the  TMOZO  mode  has  opposite  polarity  in  the  outer  and  the 

inner  radius,  this  mode  decelerates  the  annular  beam  and  accelerates  the 
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pencil  beam  if  the  beams  enter  the  cavity  in  the  same  phase.  As  the  driver 
beam  loses  energy  and  the  pencil  beam  gains  energy,  however,  the  two 
beams  will  drift  out  of  phase  with  one  another,  and  acceleration  will 
cease.  This  phase  slippage  has  a  simple  cure  [1]:  At  the  locations  where 
the  bunches  of  both  beams  enter  the  cavity  with  the  same  phase,  place  the 
annular  beam  radius  (ro)  outside  the  the  field  null  (a)  of  the  TMqzO  mode. 

When  the  bunches  of  both  beams  arrive  at  the  cavity  180°  out  of  phase, 
place  ro  inside  a.  The  radius  modulation  can  be  accomplished  by  varying 

the  external  solenoidal  magnetic  field.  Figure  2  shows  an  example  of 
acceleration  of  a  test  electron  on  the  secondary  beam  from  700  keV  to  5 
MeV  in  approximately  one  meter  of  accelerating  structure;  the  driver  beam 
is  available  from  MELBA  [3]. 

The  energy  gain  per  cavity  by  the  secondary  beam  in  a  twobetron  is 

<  Es  >=  (16.3keV)xQx(I(i  /  IkA)  x(A/  a).  (1) 

where  Q  is  the  cavity  quality  factor  of  the  TMqzO  mode  and  A  is  the 
amplitude  in  the  modulation  of  the  primary  beam  radius.  The  transformer 
ratio,  R,  which  is  the  ratio  of  the  energy  gain  in  the  secondary  beam  to  the 
energy  loss  in  the  primary  beam,  is 


R  =0.803(a/A).  (2) 

The  maximum  amount  of  secondary  beam  current.  Is,  that  can  be 
accelerated  is  limited 
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Beam  Loading 

With  radius  modulation,  the  driver  beam  must  necessarily  pass 
through  the  field  null  of  the  TMo20  mode.  When  this  occurs,  the  driver 

beam  cannot  perform  work  on  the  mode  (J*E=0),  whereas  the  secondary 
beam  does  work  against  the  mode,  causing  the  secondary  beam  to 
decelerate.  By  relaxing  the  assumption  of  a  test  particle  on  the  secondary 
beam,  we  study  the  effect  of  beam-loading.  The  result  is  shown  in  Fig.  3 
where  we  see  that,  as  Is  is  increased,  the  final  beam  energy  in  90  cm  is 

reduced. 


Beam  Breakup  Instability 

Another  theoretical  issue  that  has  been  addressed  is  the  primary 
beam  instabilities.  A  beam  with  sizable  current  traveling  down  a  series 
of  pill  box  cavities  is  vulnerable  to  a  host  of  instabilities.  Of  most 
concern  is  the  beam  breakup  instability  (BBU).  However,  we  discovered  [4] 
that  annular  beams  may  be  up  to  6  times  more  stable  than  pencil  beams 
against  the  deflecting  dipole  TMno  mode.  That  is,  the  annular  beam  can 

carry  6  times  as  much  current  as  a  pencil  beam  and  suffer  the  same  BBU 
growth.  For  the  parameters  of  the  MELBA  example,  the  number  of  e-folds 
for  BBU  growth  is  1.8  for  a  500  ns  beam  [Fig.  2]. 

Other  beam  instabilities  such  as  the  Robinson  instability  and 
Wakefield  effects  do  not  pose  a  threat  to  the  Twobetron  concept  [5]. 
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Simulation 


Several  primary  beam  issues  have  been  addressed  with  simulation. 

In  the  simulation,  40  cavities  of  the  accelerating  structure  have  been 
simulated  with  the  MAGIC  code,  which  is  a  2.5  dimensions,  fully 
relativistic,  electromagnetic  particle-in-cell  code  [6].  In  all  cases,  the 
MELBA  parameters  [Fig.  2]  were  used.  The  first  case  involved  the  beam 
traveling  down  the  accelerating  structure  to  test  that  the  modulation  is 
preserved,  and  that  no  virtual  cathodes  were  formed.  As  can  be  seen  from 
Fig.  4,  the  beam  travels  the  length  of  the  structure  with  no  disruption  and 
the  current  modulation  remains  intact.  If  we  look  at  the  cavity  fields, 
however,  we  see  that  the  structure  is  experiencing  serious  mode 
competion  with  both  the  TMqi  0  and  TMqsO  competing  with  the  desired 
TM020  mode  [Fig.  5].  This  agrees  with  intuition  since  the  beam  is 
placed  near  the  field  null  of  the  TMo20  mode,  which  approximately 
corresponds  to  the  maximum  of  the  TMOSO  mode.  Clearly,  mode  competion 
is  a  challenging  issue  to  the  Twobetron  concept.  Furthermore,  the  phase 
space  plots  of  the  beam's  energy  [Fig.  6]  show  that  the  driver  beam  is  both 
accelerated  and  decelerated  by  the  structure.  Further  tests  show  that  the 
twobetron,  in  its  present  configuration  is  a  traveling  wave  device.  This 
implies  that  a  traveling  wave  formulation  must  be  developed  for  the 
Twobetron,  or  work  must  be  done  to  isolate  the  cavities 
electromagnetically  [5]. 
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Driver  Beam  Modulation 

Simulation  was  also  performed  on  primary  beam  modulation  for  the 
MELBA  example.  The  geometry  [Fig.  7]  is  similar  to  the  existing  methods 
of  intense  beam  modulation  [7]  used  in  relativistic  klystron  amplifiers 
[RKA's].  For  a  low  current  driver  beam  (500A),  the  modulation  was  on  the 
order  of  30%  [Fig.  8].  This  was  found  to  agree  with  theory  as  there  is 
simply  not  enough  space  charge  in  the  beam  for  full  modulation  at 
convenient  levels  of  external  rf  drive.  We  therefore  decided  to  study  the 
effect  of  adding  a  feedback  loop  to  increase  beam  modulation  [Fig.  9].  This 
concept  is  similar  to  experiments  performed  recently  at  Phillips 
Laboratory  [8]  on  injection-locked  Relativistic  Klystron  Oscillators  (RKO), 
except  that  the  Phillips  Lab  experiment  utilized  feedback  from  the 
cavities  being  close  to  each  other  for  effective  electromagnetic  coupling. 

We  analyzed  the  above-mentioned  injection-locked  RKO  with  a 
lumped  circuit  model  of  the  cavities  [9],  including  the  finite  travel  time 
for  the  beam  between  the  cavities  and  the  phase  behavior  of  the  fields  in 
the  cutoff  region.  Even  though  the  model  is  fairly  simple,  it  showed 
excellent  agreement  with  experimental  results[8].  The  major  conclusion 
is  that  coupling  could  destabilize  the  system  to  provide  current 
modulation  in  a  far  shorter  drift  space  (13  cm  vs.  30  cm)  than  a 
conventional  RKA  with  a  lower  gap  voltage  on  the  booster  cavity,  thus 
avoiding  breakdown  and  virtual  cathode.  The  drawback,  when  applied  to 
the  twobetron  concept,  is  that  there  is  a  current  threshold  for  the 
operation  In  the  injection-locked  mode,  which  for  the  parameters  studied 
was  7.2kA.  Such  a  high  current  may  make  the  primary  beam  vulnerable  to 
BBU  if  the  accelerating  structure  is  excessively  long. 
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Experiments  and  Simulations 
Beam  Transport  simulations 

The  SLAC  Electron  Trajectory  Program[10]  (EGUN)was  used  in 
beam  transport  studies  to  determine  the  alpha  spread  of  an  annular 
beam  as  well  as  beam  transport  characteristics.  For  the  first  set  of 
measurements  the  4.5  cm  radius  cathode  was  modeled.  The 
parameters  used  for  EGUN  simulations  were  as  follows;  cathode 
voltage  750  keV,  diode  radius  20  cm,  drift  tube  radius  7.5  cm,  A-K 
gap  8.2  cm.  The  cathode  was  modeled  to  emit  between  4.25  and  4.75 
cm.  This  experiment  was  run  using  glyptal  enamel  to  prevent 
emission  from  other  surfaces.  In  the  computer  simulation  the  first 
2  cases  used  a  magnetic  field  similar  to  that  measured  in  the 
experiment.  Case  1  results  are  shown  in  figure  1 0.  These  show  a 
transported  current  of  3.2  kA  with  4.07  kA  emitted  from  the 
cathode.  As  figures  1 0  and  11  show,  there  is  some  beam  scraping 
going  through  the  anode  aperture.  The  average  beam  alpha 
(Vperp/Vparallel)  was  determined  to  be  approximately  0.45.  Case  2 
shown  in  figure  1 1  shows  the  same  model  with  a  grounded  plane 
across  the  anode  surface.  This  plane  simulates  aluminized  mylar 
used  in  the  experiment.  With  the  grounded  plane,  only  2.28kA  of 
current  is  seen  in  the  drift  tube  as  opposed  to  3.2  kA  of  case  1 .  But 
as  can  be  seen  from  EGUN  simulations  the  beam  is  much  more 
focused  in  the  second  case.  The  beam  needs  to  be  focused  to  pass 
through  accelerating  cavities  with  slots  that  are  0.5  cm  wide.  The 
average  beam  alpha  for  the  second  cases  was  0.45,  the  same  as  case 
1 .  Experimentally  currents  ranging  on  the  order  of  1 .8kA  were 
measured  using  the  configuration  of  case  2. 
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A  more  realistic  model  of  a  transport  experiment  is  shown  in 
figure  12.  This  model  used  a  cathode  that  emitted  at  3.5  +/-  0.25 
cm.  The  cathode  emitted  into  a  4  inch  drift  tube  after  passing 
through  the  diode.  The  A-K  gap  was  set  to  1 1  cm.  The  Magnetic  field 
of  the  diode  was  approximately  IkGauss  and  the  drift  tube  raises  the 
magnetic  field  to  2  kGauss.  The  diode  current  is  2.8kA  which  is 
fully  transported  down  the  drift  tube  in  the  simulation.  The  average 
beam  alpha  measured  was  0.6.  In  experimental  comparison  2kA  was 
transported  through  the  drift  tube,  but  the  initial  diode  current  was 
on  the  order  of  5-1  OkA. 

One  final  case  shown  in  figure  1 3  shows  the  same  parameters 
as  in  figure  12,  but  the  cathode  is  now  modeled  with  a  radius  of  2.25 
cm  with  an  emitting  surface  of  2.0  to  2.5  cm.  The  emitting  surface 
is  also  flat  vs  curved  as  in  the  other  cathodes.  The  transported  tube 
current  produced  in  EGUN  showed  1 .84kA  transported.  The  average 
beam  alpha  in  the  simulation  was  0.2.  Experimentally,  the  tube 
current  at  2kGauss  has  been  measured  to  range  between  1  and  3  kA, 
but  diode  current  ranges  between  6  and  1 0  kA. 

Experimental  Results: 

Accelerating  Cavity  Design  and  Cold  Tests 

Using  parameters  listed  in  [1],  accelerating  cavities  were 
designed.  Since  the  RF  electric  field  null  associated  with  an  axial 
electric  field  of  the  TMqzO  cavity  mode  occurred  at  3.146  cm  [1], 
cavity  slots  were  placed  at  3.35  +/-  0.25  cm.  The  cavity  radius  was 
7.24cm.  The  cavity  was  1  cm  thick.  Dimensions  of  the  accelerating 
cavity  are  shown  in  figure  14.  A  different  cavity  without  slots  was 
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also  constructed.  Cold  tests  of  the  cavity  with  slots  resulted  in 
frequency  measurements  of  figure  1 5.  The  cold  test  utilized  an  HP 
8510  Network  analyzer.  Port  one  of  the  two  port  network  was 
connected  to  a  loop  antenna  placed  inside  one  of  the  slots.  A 
monopole  antenna  was  placed  in  the  center  of  the  cavity.  The  first 
peak  in  Figure  1 5  occurs  at  3.47GHz.  This  result  can  be  accounted 
for  by  one  of  two  modes,  either  the  TMzi  (3.39  GHz)  or  the  TM02 
(3.64  GHz).  The  second  peak  occurs  at  4.0  GHz;  this  is  due  to  the 
TM3i(4.21  GHz).  And  the  third  peak  occurs  at  4.39GHz  which  is  due 
to  the  TMi  2  mode.  The  TM02  mode  had  a  Q  of  1 1 .  The  TM31  mode 
appeared  dominant  with  a  Q  of  80.  This  could  be  due  to  the  design 
with  three  slots. 

The  cavity  without  slots  for  a  beam  was  built  with  two  loop 
antennas  placed  at  the  radius  of  the  beam  slots  (r=3.15  cm).  Port 
one  of  the  network  analyzer  was  connected  through  a  power  divider 
to  the  two  antennas.  A  monopole  antenna  was  placed  in  the  center  of 
the  cavity  and  connected  to  port  two  of  the  network  analyzer.  Figure 
1 6  shows  the  results.  The  TM02  mode  demonstrated  a  weak 
response  and  occurred  at  3.76GHz,  but  the  dominant  response  was 
due  to  the  TM4i  (5.004GHz)  peak  which  occurred  at  4.88  GHz. 
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Annular  e-Beam  Transport  Experiments: 


Transport  experiments  have  been  conducted  using  three 
different  cathodes.  The  cathode  radii  are  2.25cm  +/-0.25  cm, 

3.55cm  +/-  0.25  cm,  and  4.5cm+/-0.25  cm.  Two  cathodes  are  shown 
in  figure  1 7.  Initial  experiments  were  conducted  with  an 
experimental  configuration  as  shown  in  figure  1 8  (without  the  brass 
cavities).  Experiments  used  a  6  inch  drift  tube  with  a  graphite 
collector  to  measure  transported  current.  A  piece  of  aluminized 
mylar  was  placed  in  the  anode  aperture  to  give  better  results  (EGUN 
simulations  showed  this  effect).  Experimental  parameters  are 
listed  in  figure  19  [16].  Using  the  cathode  with  a  radius  of  4.5cm 
and  an  empty  drift  tube,  measurements  were  made;  this  was 
followed  by  placing  a  brass  plate  with  three  slots  (radius  2.15cm+/- 
0.25  cm)  in  the  drift  tube  between  the  anode  and  the  graphite  plate. 
Data  recorded  are  shown  in  figures  20  and  21.  The  brass  plate  was 
replaced  by  a  pillbox  cavity  with  three  slots  (radius  2.25cm+/-  0.25 
cm)  for  measurements  recorded  in  figure  22.  Two  pillbox  cavities 
were  use  to  measure  transported  current  in  figure  23.  Without 
placing  anything  in  the  drift  tube,  the  maximum  transported  current 
was  approximately  1.8kA.  With  a  brass  plate,  the  maximum 
transported  current  was  900A  at  a  magnetic  field  of  SkGauss. 

Placing  a  pillbox  cavity  in  the  drift  tube  reduced  the  maximum 
transported  current  to  680A  at  3.7kGauss.  With  two  pillbox 
cavities,  the  maximum  transported  current  was  less  than  300A  at 
3.9kGauss.  Figure  24  shows  the  results  of  one  MELBA  shot  with  a 
pillbox  cavity  in  the  drift  tube. 
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The  cathode  with  a  3.55cni  radius  and  the  cathode  with  the 
2.25  cm  radius  were  used  with  a  4  inch  drift  tube  for  transport 
current  measurements.  Figure  25  shows  the  results  utilizing  the 
3.55  cm  radius  cathode.  These  measurements  did  not  use  aluminized 
mylar.  Figure  26  shows  the  transported  current  using  the  2.25  cm 
radius  cathode. 

Figure  27  shows  the  final,  proposed  Twobetron  experimental 
design  showing  the  orientation  of  the  accelerating  cavities  with  the 
primary  and  secondary  beams. 


Conclusions 

The  Twobetron  concept  is  a  viable,  intriguing  candidate  for 
industrial  and  medical  accelerators.  While  challenges  exist,  the  potential 
to  upgrade  pulsed  power  systems  and  develop  new  compact  accelerators 
argues  for  the  continued  development  of  this  two-beam  accelerator 
scheme. 
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A  new  configuration  is  proposed  wherein  a  low-current  beam  is  accelerated  to  high 
amps,  tens  of  McV)  by  a  driver  beam  of  high  current  and  low  energy  (a  few  kiloam^,  <  MeV).  The 
annular  driver  beam  excites  the  TMom  cavity  mode  of  an  accelerating  structure  which  transfers  its  rf 
power  to  the  on-axis  secondary  beam.  Systematic  variation  of  the  driver  beam  radius  provides  the  secon¬ 
dary  beam  with  phase  focusing  and  adjustable  acceleration  gradient.  A  proof-of-principle  experiment  is 

suggested. 
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Compact  electron  and  ion  accelerators  in  the  10  MeV 
range  have  a  wide  range  of  applications,  such  as  treat¬ 
ment  of  bulk  materials,  activation  analysis,  and  medical 
radiation  sources.  To  achieve  such  an  energy  at  moderate 
levels  of  current  (tens  of  amps)  requires  considerable 
power,  and  a  natural  candidate  for  a  driver  is  the  pulse 
power  system  [1,2].  Intense  annular  electron  beams  (a 
few  kiloamps,  <  1  MeV)  extracted  from  such  a  system 
have  been  modulated  efficiently,  and  the  current  modula¬ 
tions  exhibit  a  high  degree  of  amplitude  and  phase  stabil¬ 
ity  [3].  These  modulated  beams  have  been  used  to  gen¬ 
erate  ultrahigh  power  microwaves  [4,5]  and  to  accelerate 
electrons  to  high  energies  [6].  They  will  be  used  as  the 
driver  in  the  two-beam  accelerator  to  be  proposed  in  this 

paper.  .  . 

Various  two-beam  accelerators  have  been  studied  in 
the  past  [6-10].  There  are  significant  differences  in  the 
present  configuration,  shown  schematically  in  Fig.  1.  The 
driver  beam  is  an  annular  beam  of  radius  ro,  carrying  an 
ac  current  /</  at  frequency  (o.  It  passes  through  an  ac¬ 
celerator  structure,  consisting  of  N  cylindrical  pillbox 
cavities.  Each  cavity  has  a  radius  6  =5.52c/<a  so  that  (o 
is  also  the  resonant  frequency  of  the  TM020  mode  of  the 
pillbox  cavity  (Fig.  I).  The  secondary  beam  is  an  on-axis 
pencil  beam,  carrying  an  ac  current  Is  (Is'^Idh  aiso  at 
frequency  ©.  Since  the  rf  electric  fields  of  the  TM020 
mode  have  opposite  signs  in  the  outer  region  and  in  the 
inner  region,  the  mode  retards  the  annular  driver  beam 
but  accelerates  the  on-axis  secondary  beam.  As  we  shall 
see,  if  the  driver  beam  radius  is  modulated  axially,  phase 
focusing  and  tunability  in  the  output  energy  of  the  secon¬ 
dary  beam  can  be  achieved.  This  is  the  crucial  feature  of 
the  present  device,  not  shared  by  the  prior  works  [6-10]. 

Thus,  without  the  use  of  rf  plumbing,  the  present 
scheme  provides  the  gradual  conversion  of  the  primary 
beam  power  to  the  secondary  beam  over  many  accelerat¬ 
ing  gaps.  Since  the  current  modulation  on  the  primary 
beam  has  been  shown  to  be  insensitive  to  the  variations  in 
the  diode  voltage  and  diode  current  [3],  the  effectiveness 
in  the  acceleration  of  the  secondary  beam  is  likewise  in- 
sensitive  to  such  variations. 

To  calculate  the  excitation  of  the  TM020  mode  by  the 


primary  beam,  and  the  resultant  acceleration  of  the 
secondary  beam  by  this  mode,  we  assume  that  the  intense 
space  charge  on  the  beam  does  not  alter  the  rf  charac¬ 
teristic  of  the  cavities  [4,1 1,12].  We  also  assume  that  the 
individual  pillbox  cavities  are  electromagnetically  isolated 
from  each  other  when  the  beams  are  absent  [13,14]. 
Since  the  cavities  are  excited  mainly  by  the  rf  current  Ij 
carried  by  the  primary  beam,  the  TM020  mode  so  excited 
always  decelerates  the  primary  beam  electrons  on  the 
average  (by  conservation  of  energy).  This  is  true  whether 
the  beam  radius  ro  is  larger  or  smaller  than  a,  where 
a  =2.405c/a)  is  the  radius  of  the  rf  electric  field  null  of 
the  TM020  mode  [Fig.  2(a)].  The  value  of  the  rf  electric 
field  at  ro  then  gives  the  deceleration  gradient.  In  terms 
of  the  relativistic  mass  factor  (/</),  the  energy  loss  by  this 
driver  beam  as  it  traverses  the  nth  cavity  is  given  by 


dn 


I  a  continuum  description.  In  Eq.  (1), 


.  ^  ^  /  w  l  /l 


(1) 

(2) 


is  the  dimensionless  parameter  that  measures  the  strength 
of  the  cavity  excitation  by  the  primary  beam, 

5-yo(<aro/c)«  -  1.249(ro-a)/a  ,  (3) 

Q  is  the  quality  factor  of  the  TM020  mode,  L  is  the  cavity 
length,  and  Jo  is  the  Bessel  function  of  the  first  kind  of 


driver  beam  (annular,  modulated)  j^((i)r/c) 


secondary  beam  (pencil,  modulated) 


FIG.  I.  Schematic  drawing  of  the  two-beam  accelerator. 
Also  shown  is  the  rf  force  profile,  Joimr/c),  associated  with  the 
axial  electric  field  of  the  TM020  cavity  mode. 
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(A) 


(B) 


FIG.  2.  (a)  Position  of  the  primary  beam  radius  ro  (ro>a) 
for  secondary  beam  acceleration  when  both  beams  enter  the 
cavity  at  the  same  phase,  (b)  Position  of  the  primary  beam  ra¬ 
dius  ro  (ro<a)  for  secondary  beam  acceleration  when  both 
beams  enter  the  cavity  at  1 80°  phase  apart. 


order  zero.  In  writing  the  last  expression  of  Eq.  (3),  we 
have  made  the  assumption  that  the  annular  beam  is  locat¬ 
ed  in  the  vicinity  of  the  rf  electric  field  null  (ro«  a). 

If  the  secondary  beam  enters  the  cavity  at  the  same 
phase  as  the  primary  beam,  the  former  will  be  accelerat¬ 
ed  if  ro  >  a,  for  in  this  case  the  rf  fields  experienced  by 
both  beams  have  opposite  polarity  [Fig.  2 (a)].  Sinw  the 
rf  electric  field  has  a  radial  dependence  of  Joiarlc),  it  is 
obvious  that  1/|5|  is  the  “transformer  ratio,”  which  is  the 
ratio  of  the  energy  gain  by  the  secondary  beam  to  the  en¬ 
ergy  loss  by  the  primary  beam,  if  both  beams  enter  the 
cavity  at  the  same  phase.  This  dependence  on  the  phase 
is  reflected  in  the  following  equation  which  describes  the 
change  in  the  relativistic  mass  factor  (y,)  of  the  secon¬ 
dary  beam  as  it  traverses  the  nth  cavity: 


-^  =  -A5cos(0,-0rf),  ^4) 

dn 

where  Os  is  the  phase  of  the  secondary  beam  bunch  and 
Od  is  the  phase  of  the  primary  beam  bunch  when  they 
enter  the  nth  cavity.  Equation  (4)  is  readily  obtained 
from  Eq.  (1)  by  noting  the  transformer  ratio  1/5  and  the 
phase  difference  mentioned  above.  Equations  (3)  and  (4) 
indeed  show  that  y,  increases  \( ro>  a  and  if  9d  —^5- 
The  secondary  beam  cannot  be  accelerated  indefinitely 
because  of  the  increase  in  the  phase  slippage  between  Od 
and  9s  downstream.  This  phase  slippage  occurs  as  the 
primary  beam  is  decelerated  and  the  secondary  beam  is 
accelerated.  Its  rate  of  increase  is  governed  by 


d(9s  ~9d)  ^  o)L 
dn  c 


Pd 


=  "^[(1  _  i/y2)  -'^2-  (1  -  i/ri)  . 

C 

(5) 


FIG  3  Evolution  of  the  relativistic  mass  factors  when  the 
driver  'beam  radius  ro  is  a  constant:  (a)  the  driver  beam  (b) 
the  secondary  beam.  Phase  slippage  prohibits  continual  ac- 
celeration  of  the  secondary  beam. 

The  effect  on  the  secondary  beam  by  this  phase  slippage 
is  illustrated  in  Fig.  3,  which  is  obtained  by  numerically 
solving  the  system  of  three  equations  1(1),  (4),  (571  in 
three  unknowns:  yj,  y„  and  9s -6,.  The  initiaUondi- 
tions  for  these  three  unknowns  are  taken  to  Iw  6,  dj  -0 
and  yrf“yj“°2.37,  corresponding  to  an  initial  energy  of 
700  keV  for  both  beams.  The  other  parameters  are 
6)/2;r“3.65  GHz,  6-7.221  cm,  L  =  1  cm,  n  =3.146  cm, 
ro= 3.322  cm,  0  =  100,  and  /rf=0.5  kA.  Since  we  have 
taken  L  =  1  cm,  the  cavity  number  n  is  also  the  axial  dis¬ 
tance  (z)  in  cm.  ,  •  •4-  i 

Figure  3(a)  shows  that  yj  decreases  from  the  initial 

value  of  2.37  to  1.24  at  «=90;  i.e.,  the  primary  beam  s 
energy  steadily  decreases  from  700  to  125  keV  after 
propagating  90  cm.  The  secondary  beam’s  energy  IFig. 
3(b)l  increases  initially,  reaching  a  maximum  value  of 
2  3  MeV  after  24  cm,  and  then  decreases  due  to  the 
phase  slippage  until  «=56,  and  oscillates  further  down¬ 
stream  as  the  phase  slippage  continues. 

The  phase  slippage  may  be  corrected  by  adjusting  the 
primary  beam’s  radius  ro-  Consider,  for  example,  the 
worst  case  of  phase  slippage  where  the  primary  beam  and 
the  secondary  beam  arrive  at  a  cavity  180“  out  of  phase, 
as  shown  in  Fig.  2(b).  If  the  primary  beam’s  radius  ro  is 
less  than  a,  it  generates  an  rf  electric  field  which  would 
retard  both  beams  during  the  time  when  the  primary 
beam  occupies  the  cavity.  However,  when  the  charge 
bunch  of  the  primary  beam  resides  in  the  cavity,  there  are 
few  particles  in  the  secondary  beam  residing  in  the  same 
cavity  because  both  beams  arrive  at  the  cavity  18()“  out 
of  phase.  By  the  time  the  charge  bunch  of  the  primary 
beam  is  about  to  leave  the  cavity,  the  rf  electric  field  is 
about  to  change  sign,  at  which  time  the  charge  bunch  of 
the  secondary  beam  is  about  to  enter  the  cavity,  whose  rf 
electric  field  then  begins  to  accelerate  the  entering  bunch 
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on  the  secondary  beam.  Thus,  the  phase  slippage  prob¬ 
lem  can  be  corrected  by  a  simple  cure:  At  the  locations 
where  the  bunches  of  both  beams  enter  the  cavity  with 
the  same  phase,  place  ro  outside  a.  When  the  bunches  of 
both  beams  arrive  at  the  cavity  1 80®  out  of  phase,  place 

ro  inside  a.  r-  a  fA\ 

Mathematically,  it  is  easy  to  see  from  Eqs.  131  and 

that  Ys  is  3  monotonically  increasing  function  of  n  if  ro  is 
tapered  in  such  a  way  that  (ro~a) cos(0s  —  0j)>O. 

The  above  idea  of  phase  slippage  correction  has  l^en 
tested  for  the  example  shown  in  Figs.  3(a)  and  3(b). 
From  that  figure,  the  phase  slippage  occurs  with  a  period 
of  the  order  of  75  cm.  Thus,  we  correct  the  primary 
beam  radius  ro  according  to 

ro(cm )  =  3. 1 46  -b  (3.322  -  3. 1 46 )  cos(2;r«/75  ) .  (6) 

Including  only  this  modification,  and  keeping  all  other 
parameters  the  same,  we  obtain  Fig.  4.  In  Fig.  4,  we  see 
that  the  primary  beam’s  energy  monotonically  decreases 
from  700  to  400  keV  over  90  cm,  whereas  the  secondary 
beam’s  energy  increases  monotonically  from  700  keV  to  a 
maximum  of  4.2  MeV  over  the  same  distance,  in  sharp 
contrast  to  Fig.  3(b).  The  loss  of  300  keV  in  the  primary 
beam  and  the  gain  of  3.5  MeV  in  the  secondary  beam 
implies  an  effective  transformer  ratio  of  about  (3.5 
MeV)/(300keV)=11.7. 

In  Fig.  4,  the  zero  slopes  in  y,  and  in  yj  occur  at  the 
axial  positions  (n)  at  which  the  driver  beam  radius  ro 
coincides  with  the  field-null  position  a.  The  slight  dip  in 
Ys  at  n  =90  only  means  that  the  primary  beam’s  radius^ro 
needs  further  adjustment  there.  If  we  write  ro~o 
-I-Acos(v'),  where  A  is  the  amplitude  and  yf  is  the  phase 
of  the  modulation  in  ro,  the  general  phase  focusing  condi¬ 
tion  reads  dyr/dn  ^diO,  -  Qd)ldn.  This  condition  is  appl¬ 
icable  when  the  two  beams  have  different  velocities.  In 
fact,  one  might  argue  that  this  technique  of  radius  modu¬ 
lation  provides  both  beams  with  self-focusing  in  phase, 
similar  to  the  self-focusing  in  synchrotrons  [151. 

The  modulation  in  the  annular  beam  radius  may  be 
readily  achieved  by  a  proper  adjustment  of  the  external 
solenoidal  magnetic  field  which  is  often  used  for  beam 
focusing  and  beam  transport  [3-6,14].  Since  the  rate  of 
change  of  energy  depends  on  the  annular  beam,xadius  ro 
[cf.  Eqs.  (I)  and  (3)1,  the  output  energy  of  the  accelerat¬ 
ed  beam  may  also  be  controlled  by  the  same  external 

magnetic  field  coils.  _ 

The  above  ideas  may  be  tested  in  a  proof-of-principle 
experiment  with  parameters  similar  to  those  used  to  pro¬ 
duce  Fig.  4.  The  primary  beam  may  be  obtained,  for  ex¬ 
ample,  from  the  Michigan  Electron  Long-Beam  Ac¬ 
celerator  (MELBA)  [16],  which  operates  with  diode  pa¬ 
rameters  of  700  keV,  current  up  to  10  kA,  and  flattop 
pulse  length  up  to  1  /is.  This  primary  beam  may  be 
modulated  using  the  proven  techniques  by  Friedman  et 
al.  [3,4,61.  Note  that  the  average  acceleration  gradient 
of  40  kV/cm  and  the  peak  acceleration  gradient  of  about 
80  kV/cm  implied  by  Fig.  4  are  well  within  the  rf  break- 

49- 


FIG.  4.  Evolution  of  the  relativistic  mass  factors  when  the 
driver  beam  radius  ro  is  modulated  to  compensate  phase  slip¬ 
page:  (a)  the  driver  beam,  (b)  the  secondary  beam. 

down  limit.  If  we  assume  an  acceleration  efficiency  of 
25%,  a  secondary  beam  of  more  than  10  A  of  current 
may  be  accelerated  to  4  MeV  in  less  than  a  meter  in  this 
proof -of-principle  experiment. 

There  are  many  issues  which  may  affect  the  eventual 
usefulness  of  the  two-beam  accelerator  concept  outlined 
above.  Chief  among  them  is  the  modification  of  the  rf 
characteristic  that  always  accompanies  an  intense  driver 
beam,  which  includes  a  detuning  of  the  structure  frequen¬ 
cy  and  a  modification  of  the  gap  transit-time  factor 
[4,1 1,12,17].  Also  of  concern  is  the  beam  breakup  insta¬ 
bility  (BBU)  on  the  driver  beam  [10,13,14,17].  However, 
we  have  recently  found  that  BBU  in  an  annular  beam 
may  be  far  less  serious  than  a  pencil  beam  [18],  and  BBU 
can  be  controlled  by  many  well-known  techniques  [19]. 
The  degree  of  coupling  among  neighboring  cavities,  espe¬ 
cially  in  the  presence  of  an  intense  beam,  remains  to  be 
studied  [20].  Although  the  driver  beam’s  radius  is  a  cru¬ 
cial  factor,  the  effects  of  the  beam’s  finite  thickness  are 
far  less  important,  according  to  our  preliminary  studies. 
We  have  also  examined  the  effects  of  the  transverse  wake 
[21]  and  of  the  longitudinal  instabilities  [22]  and  found 
that  they  are  not  serious,  at  least  for  the  parameters  used 
in  the  above  numerical  example,  assuming  a  solenoidal 
field  of  10  kG  in  the  accelerating  structure. 

In  summary,  we  propose  a  novel  scheme  which  has  the 
potential  of  converting  many  existing  pulse  power  systems 
into  compact  rf  accelerators  that  are  suitable  for  industri¬ 
al  and  medical  applications.  The  driver  beam  is  a  modu¬ 
lated  intense  relativistic  electron  beam  of  annular  shape 
and  low  energy  (  <  1  MeV).  The  secondary  beam  is  an 
on-axis  pencil  beam.  The  secondary  beam  may  reach  an 
energy  up  to  10  MeV  in  1  to  2  m.  Phase  focusing  and  en¬ 
ergy  tunability  of  the  accelerated  beam  may  be  provided 
by  an  external  magnetic  field,  which  controls  the  radius 
of  the  primary  beam.  While  we  have  in  this  paper  con- 
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centrated  only  on  electron  acceleration  in  the  10  MeV 
range,  it  is  intriguing  to  speculate  on  the  potential  of  us¬ 
ing  this  technique  (a)  to  accelerate  ions  to  tens  of  MeV, 
and  (b)  to  accelerate  electrons  to  ultrahigh  energy  using 
superconducting  cavities  Icf.  Eq.  (2)1  and  higher  energy 
driver  beams. 

We  thank  John  W.  Luginsland  for  assistance  in  the 
preparation  of  this  manuscript.  This  work  was  supported 
by  SDIO-BMD/IST/ONR. 
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Abstract.  A  new  configuration  is  analyzed  wherein  a  low  current  beam  is  accelerated  to  high 
energies  (lO's  of  amps,  lO's  of  MeV)  by  a  driver  beam  of  high  current  and  low  energy  (a  few 
kiloamps,  <  1  MeV).  The  annular  driver  beam  excites  the  TMo20  cavity  mode  of  an  acceleraUng 
structure  which  transfers  its  rf  power  to  the  on-axis  secondary  beam.  Systematic  variation  of  the 
driver  beam  radius  provides  the  secondary  beam  with  phase  focusing  and  adjustable  acceleration 
gradient.  A  proof-of-principle  experiment  is  suggested.  Various  issues,  such  as  the  scaling  laws, 
transverse  and  longitudinal  instabilities,  rf  coupling  among  cavities,  etc.,  are  examined. 


INTRODUCTION 


Two-beam  accelerators  have  been  studied  extensively  in  the  high  energy 
physics  community  [1].  This  paper  concentrates  on  the  10  MeV  range.  Compact 
electron  and  ion  accelerators  in  this  energy  range  have  a  wide  range  of 
applications,  such  as  treatment  of  bulk  materials,  activation  analysis,  and  medical 
radiation  sources.  To  achieve  such  an  energy  at  moderate  levels  of  current  (tens 
of  amps)  requires  considerable  power,  and  a  natural  candidate  for  a  driver  is  the 
pulse  power  system.  Intense  annular  electron  beams  (multi-kiloamps,  <  1  MeV) 
extracted  from  such  a  system  have  been  modulated  efficiently,  and  the  current 
modulations  exhibit  a  high  degree  of  amplitude  and  phase  stability  [2].  Their 
successful  applications  [3]  in  ultra-high  power  microwave  generation  and  in 
particle  acceleration  have  motivated  us  to  use  them  as  drivers  in  a  novel  two-beam 
accelerator  [4],  termed  "twobetron"  hereafter. 

The  twobetron  is  shown  schematically  in  Fig.  1.  The  driver  beam  is  an 

annular  beam  of  radius  ro,  carrying  an  AC  current  Id  at  frequency  (O.  It  passes 
through  an  accelerator  structure,  consisting  of  N  cylindrical  pillbox  cavities.  Each 
cavity  has  a  radius  b  =  5.52  c/m  so  that  ©  is  also  the  resonant  frequency  of  the 
TMo20  mode  of  the  pillbox  cavity  [Fig.  1].  The  secondary  beam  is  an  on-axis 
pencil  beam,  carrying  an  AC  current  Is  (Is  « Id),  also  at  frequency  ©.  Since  the 
rf  electric  fields  of  the  TM020  mode  have  opposite  signs  in  the  outer  region  and  in 
the  inner  region,  the  mode  retards  the  annular  driver  beam  but  accelerates  the  on- 
axis  secondary  beam.  If  the  driver  beam  radius  is  modulated  axially ,  phase 
focusing  and  tunability  in  the  output  energy  of  the  secondary  beam  can  be 
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achieved  This  phase  focusing  technique  thus  also  offers  the  possibility  of  using 
much  lower  driver  beam  voltage  (e.g.,  100  keV),  a  distinct  advantage  m  many 

applications. 


driver  beam  (annular,  modulated)  J^(cor/c) 


0 


To 


TM020  Mode 


secondary  beam  (pencil,  modulated) 


Figure  1.  Schematic  drawing  of  the  two-beam  accelerator.  Also  shown  is  the  rf 
force  profile,  Jo(CDr/c),  associated  with  the  axial  electric  field  of  the  TM020  cavity 
mode. 


ACCELERATION  MECHANISM 

To  calculate  the  excitation  of  the  TM020  mode  by  the  primary  beam,  and  the 
resultant  acceleration  of  the  secondary  beam  by  this  mode,  we  assume  that  the 
intense  space  charge  on  the  beam  does  not  alter  the  rf  characteristic  of  the 
cavities.  We  also  assume,  for  the  time  being,  that  the  individual  pillbox  cavities 
are  electromagnetically  isolated  from  each  other  when  the  beams  are  absent. 
Since  the  cavities  are  excited  mainly  by  the  rf  current  Id  carried  by  the  primary 
beam,  the  TM020  mode  so  excited  always  decelerates  the  primary  beam  electrons 
on  the  average  (by  conservation  of  energy).  This  is  true  whether  the  beam  radius 
ro  is  larger  or  smaller  than  a,  where  a  =  2.405  c/O)  is  the  radius  of  the  rf  electric 
field  null  of  the  TM020  mode  [Fig.  2a].  The  value  of  the  rf  electric  field  at  ro  then 
gives  the  deceleration  gradient.  In  terms  of  the  relativistic  mass  factor  (y^),  the 
energy  loss  by  this  driver  beam  as  it  traverses  the  n-th  cavity  is  given  by 
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(1) 


dn 


in  a  continuum  description.  In  Eq.  (1), 

A  =  0.066(a)L  /  c)Q(I^j  /  IkA)  (2) 

is  the  dimensionless  parameter  that  measures  the  strength  of  the  cavity  excitation 
by  the  primary  beam, 

5  =  J^(o)r^ /c)==-1.249(r^ -a)/a,  (3) 

Q  is  the  quality  factor  of  the  TM020  mode,  L  is  the  cavity  length,  Jq  is  the  Bessel 
function  of  the  first  kind  of  order  zero.  In  writing  the  last  expression  of  Eq.  (3), 
we  have  made  the  assumption  that  the  annular  beam  is  located  at  the  vicinity  of 
the  rf  electric  field-null  (  ro  ~  s)- 
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Figure  2. 

(a)  Left.  Position  of  the  primary  beam  radius  ro  (  Tq  >  a)  for  secondary  beam 

acceleration  when  both  beams  enter  the  cavity  at  the  same  phase. 

(b)  Right.  Position  of  the  primary  beam  radius  ro  (  ro  <  a)  for  secondary  beam 

acceleration  when  both  beams  enter  the  cavity  at  180  apart  phase. 
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If  the  secondary  beam  enters  the  cavity  at  the  same  phase  as  the 
beam,  the  former  will  be  accelerated  if  Tq  >  a,  for  m  this  case  the  rf  fields 
experienced  by  both  beams  have  opposite  polarity  [Fig.  2a].  Since  the  rf  electnc 
field  has  a  radial  dependence  of  Jo(®r/c),  it  is  obvious  that  1/  |5|  is  the 
"transformer  ratio",  which  is  the  ratio  of  the  energy  gain  by  the  secondary  beam  to 
the  energy  loss  by  the  primary  beam,  if  both  beams  enter  the  cavity  at  the  s^ 
phase.  This  dependence  on  the  phase  is  reflected  in  the  following  equation  which 
describes  the  change  in  the  relativistic  mass  factor  (Ys)  of  secondary  beam  as 
it  traverses  the  n-th  cavity: 

^  =  -A5-cos(es-ed), 
dn 

where  0s  is  the  phase  of  the  secondary  beam  bunch  and  Gd  is  the  phase  of  the 
primary  beam  bunch  when  they  enter  the  n-th  cavity.  Equation  (4)  is  readily 
obtained  from  Eq.  (1)  by  noting  the  transformer  ratio  1/5  and  the  phase  difference 
mentioned  above.  Equations  (3)  and  (4)  indeed  show  that  Ys  increases  if  ro  >  a 

and  if  6^  =  Gg. 

The  secondary  beam  cannot  be  accelerated  indefinitely  because  of  the 
increase  in  the  phase  slippage  between  0^  and  Gg  downstream.  This  phase 
slippage  occurs  as  the  primary  beam  is  decelerated  and  the  secondary  beam  is 
accelerated.  Its  rate  of  increase  is  governed  by 

d(0s  -Qd) 

dn  c  l^pg  J  c 

The  effect  on  the  secondary  beam  by  this  phase  slippage  is  illustrated  in  Fig.  3, 
which  is  obtained  by  numerically  solving  the  system  of  three  equations  [(1),  (4), 
(5)]  in  three  unknowns:  Yd’  Ts’  ^s  “  ^d'  ^he  initial  conditions  for  these  three 
unknowns  are  taken  to  be:  Gg  -  0^  =  0  ,  Yd  =  Ts  =  2-  37,  corresponding  to  an 

initial  energy  of  700  keV  for  both  beams.  The  other  parameters  are  :  (o/27t  =  3.65 
GHz,  b  =  7.221  cm,  L  =  1  cm,  a  =  3.146  cm,  ro  =  3.322  cm,  Q  =  100,  Id  =  0.5  kA. 
Since  we  have  taken  L  =  1  cm,  the  cavity  number  n  is  also  the  axial  distance  (z)  in 
cm. 

Figure  3a  shows  that  decreases  from  the  initial  value  of  2.37  to  1.24  at 

n=90,  i.e.,  the  primary  beam's  energy  steadily  decreases  from  700  keV  to  125  keV 
after  propagating  90  cm.  The  secondary  beam's  energy  [Fig.  3b]  increases 
initially,  reaching  a  maximum  value  of  2.3  MeV  after  24  cm,  and  then  decreases 
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due  to  the  phase  slippage  until  n=  56,  and  oscillates  further  downstream  as  the 
phase  slippage  continues. 


Figure  3.  Evolution  of  the  relativistic  mass  factors  when  the  driver  beam  radius 
ro  is  a  constant:  (a)  the  driver  beam,  (b)  the  secondary  beam.  Phase  slippage 
prohibits  continual  acceleration  of  the  secondary  beam. 
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The  phase  slippage  may  be  corrected  [4]  by  adjusting  the  primary  beam's 
radius  tq.  Consider,  for  example,  the  worst  case  of  phase  slippage  where  the 

primary  beam  and  the  secondary  beam  arrive  at  a  cavity  180®  out  of  phase,  as 
shown  in  Fig.  2b.  If  the  primary  beam's  radius  ro  is  less  than  a,  it  generates  an  it 
electric  field  which  would  retard  both  beams  during  the  time  when  the  primary 
beam  occupies  the  cavity.  However,  when  the  charge  bunch  of  the  primary  be^ 
resides  in  the  cavity,  there  are  few  particles  in  the  secondary  beam  residing  m  the 

same  cavity  because  both  beams  arrive  at  the  cavity  180®  out  of  phase.  By  the 
time  the  charge  bunch  of  the  primary  beam  is  about  to  leave  the  cavity,  the  rf 
electric  field  is  about  to  change  sign,  at  which  time  the  charge  bunch  of  the 
secondary  beam  is  about  to  enter  the  cavity,  whose  rf  electric  field  then  bepns  to 
accelerate  the  entering  bunch  on  the  secondary  beam.  Thus,  the  phase  slipage 
problem  can  be  corrected  by  a  simple  cure:  At  the  locations  where  the  bunches  of 
both  beams  enter  the  cavity  with  the  same  phase,  place  ro  outside  a.  When  the 

bunches  of  both  beams  arrive  at  the  cavity  180®  out  of  phase,  place  ro  mside  a. 

Mathematically,  it  is  easy  to  see  from  Eqs.  (3)  and  (4)  that  Yg  ^ 
monotonically  increasing  function  of  n  if  ro  is  tapered  in  such  a  way  that 
(ro  -  ajcosCOg  -  O^j)  >  0.  Thus,  if  we  write  the  primary  beam  radius  ro  as 

r^ -a  =  Acos(\j/), 

where  A  is  the  amplitude  and  X)/  is  the  phase  of  the  modulation  in  ro,  the  general 
phase  focusing  condition  reads 

d\|/  /  dn  =  d(0s  -  6(j )  /  dn. 


Including  this  modification,  with  A=  3.322  cm  -  3.146  cm  =  0.176  cm,  and 
keeping  all  other  parameters  the  same  as  in  Fig.  3,  we  obtain  Fig.  4.  In  Fig.  4,  we 
see  that  the  primary  beam's  energy  monotonically  decreases  from  700  keV  to  415 
keV  over  90  cm,  whereas  the  secondary  beam's  energy  increases  monotonically 
from  700  keV  to  a  maximum  of  4.8  MeV  over  the  same  distance,  in  sharp  contrast 
to  Fig.  3b.  The  loss  of  285  keV  in  the  primary  beam  and  the  gain  of  4.1  MeV  in 
the  secondary  beam  implies  an  effective  transformer  ratio  of  about  4.1  MeV/285 

keV=14.4. 

The  modulation  in  the  annular  beam  radius  may  be  readily  achieved  by  a 
proper  adjustment  of  the  external  solenoidal  magnetic  field  which  is  often  used 
for  beam  focusing  and  beam  transport.  Since  the  rate  of  change  of  energy 
depends  on  the  annular  beam  radius  ro  [cf.  Eqs.  (1),  (3)],  the  output  energy  of  the 
accelerated  beam  may  also  be  controlled  by  the  same  external  magnetic  field 

coils.  .  .  .  .  , 

The  above  ideas  may  be  tested  in  a  proof-of-principle  experiment  with 

parameters  similar  to  those  used  to  produce  Fig.  4.  The  primary  beam  may  be 


oWained,  tor  example,  from  the  Michigan  Electron  Long.BeamA^^^^ 
(MELBA  (51),  which  operates  with  diode  parameter  of  700  keV,  current  up  to  1 
kA.  and  flattop  pulse  length  up  to  1  (ts.  This  primaty  beam  may  be  modulated 
usiug  the  proven  techniques  by  Friedman  et  al.  (2],  Note  that  the  avemge 
acceleration  gradient  of  45  kV/cm  and  the  peak  accelerahon  gradient  of  about  80 
kV/cm  implied  by  Fig.  4  are  weU  within  the  rf  breakdown  hmiL  If  we  ^sime  m 
acceleration  efficiency  of  25  per  cent,  a  secondary  beam  of  more  than  8  amps  of 
current  may  be  accelerated  to  5  MeV  in  less  than  a  meter  m  this  proof-of-pnnciple 

experiment. 


Figure  4.  Evolution  of  the  relativistic  mass  factors  when  the  driver  beam  radius 
ro  is  modulated  to  compensate  phase  slippage;  (a)  the  driver  beam,  (b)  the 

secondary  beam. 
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ISSUES 


We  shall  address  some  of  the  issues  which  may  affect  the  eventual 
usefulness  of  the  twobetron.  They  concern  the  scaling  laws,  beam  breakup 
instabilities,  wakefield  effects,  effects  of  finite  beam  thickness,  coupling  among 
cavities,  modification  of  rf  characteristic  by  the  intense  driver  beam,  and  the 
integrity  of  the  primary  beam  modulation  in  the  accelerating  structure. 

Scaling 

The  average  energy  gain  per  cavity  by  the  secondary  beam  in  the  twobetron 
is 

<  Es  >=  (16.3keV)  x  Q  x  (I^  /  IkA)  x  (A  /  a) .  (^) 

The  transformer  ratio,  R,  which  is  the  ratio  of  the  energy  gain  in  the  secondary 
beam  to  the  energy  loss  in  the  primary  beam,  is 

R  =  0.803(a/A). 

The  maximum  amount  of  secondary  beam  current.  Is,  that  can  be  accelerated  is 
limited  to 

^  2R 

Given  a  driver  beam  current,  we  cannot  make  the  acceleration  gradient  (i.e., 
<Es>)  excessively  high  by  using  a  very  high  Q  cavity.  A  practical  limit  on  <Es> 

is  set  by  rf  breakdown  in  the  cavities  [6].  .  , .  . 

In  general,  a  transformer  ratio  R  of  the  order  of  10  seems  achievable.  A  two- 
stage  twobetron,  in  which  the  accelerated  beam  of  the  first  stage  is  used  as  the 
driver  beam  in  the  second  stage,  will  provide  voltage  multiplication  by  a  factor  of 
100,  while  the  output  current  is  correspondingly  much  reduced  [cf.  Eq.  10]. 

Primary  Beam  Instabilties 

The  intense  driver  beam  passing  through  a  sequence  of  cavities  is  highly 
vulnerable  to  beam  breakup  instabilities  (BBU  [7-10]).  However,  we  have 
recently  found  that  BBU  in  an  annular  beam  may  be  far  less  serious  than  a  pencil 
beam  [8].  Specifically,  in  the  absence  of  other  stabilizing  mechanisms  such  as 
stagger  tune  and  betatron  frequency  spread,  the  number  of  e-folds,  (A/,  in  BBU 
amplitude  growth)  over  an  accelerator  of  length  z  is  given  by  [9] 
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iPdC^OpJ 

where  ©b  is  the  frequency  of  the  deflecting  mode  with  quality  factor  Qb,  cop  is  the 
betatron  frequency  associated  with  the  focusing  field,  £  =  0.0041(pd/Yd)(Id/lkA)  is 
the  dimensionless  coupling  constant,  and  pd.  Yd  are  defined  in  Eq.  5.  For  a 
solenoidal  magnetic  field  of  10  kG  and  a  dipole  mode  Qb  of  100,  N  -  1.8  for  a 
500  ns,  0.5  kA  driver  beam  in  a  90  cm  accelerator  structure,  as  in  the  numerical 

conjecture  that  the  longitudinal  (Robinson-like)  instability  [10]  probably 
is  not  important  for  the  twobetron,  at  least  in  the  proposed  proof-of-princip  e 
experiment  Unlike  a  circular  accelerator,  the  present  scheme  is  single-pass.  Its 
acceleration  length  is  quite  short,  its  length  is  only  slightly  over  one  wavelength  in 
the  radius  modulation.  Moreover,  the  dnve  frequency  may  be  adjusted  to  be  on 
the  "right  side"  of  the  structure  frequency  to  avoid  the  Robinson-like  instability. 

In  a  preliminary  particle  simulation,  we  find  that  the  current  modulation  is 
preserved  on  the  primary  beam,  after  it  is  made  to  propagate  through  the 
federating  structure,  using  the  beam  and  structure  parameters  that  are  being 
planned  for  a  proof-of-principle  experiment. 

Wakefield  Effects 

We  have  also  examined  the  effects  of  the  transverse  wake  on  the  driver 
beam,  and  found  that  a  nominal  value  of  solenoidal  field  Bo  -  10  kG  would 
render  the  effects  of  the  transverse  wakefield  on  the  driver  beam  unimportant. 
Specifically,  under  the  condition  n»®/Yd(l  +  Pd)  , where  Q  is  the 
nonrelativistic  cyclotron  frequency  associated  with  Bo  and  the  other  symbols  are 
the  same  as  in  Eq.  (5),  the  electron  motion  is  adiabatic  along  the  composite  (DC  + 
rf)  magnetic  field  line.  The  maximum  angular  displacement,  from  the  mean,  is 

estimated  to  be 

4  =  0. 52(c  /  ©)(Ea  /  cBo  )Pd  /  (1  -  Pd )  ^  ^ 

where  Ea  is  the  maximum  accelerating  electric  field  experienced  by  the  secondary 
beam.  The  maximum  radial  displacement  is 

£r=^Q^L^^m->  ^  ^ 

where  =271^^1 is  the  axial  wavelength  associated  with  the 

modulation  in  the  driver  beam  radius.  For  the  parameters  used  in  the  numerical 
example,  £q  <  0.2cm,  and  ^  0.0058cm.  The  spread  in  momentum,  dp,  in  the 


b  ^ 
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driver  beam  may  introduce  a  variation  in  its  annular  beam  radius,  dro-  It  is 
estimated  that 


dro  <  [2(4.  /  4n)^  +  0.083>.LEa  '  cBoJdp  /  p  , 


(14) 


where  A  is  the  amplitude  of  the  modulation  in  the  driver  bearn  radius.  Using  the 
parameters  in  the  numerical  example  [Fig.  4],  we  find  dro  <  0.0061  cm  if  dp/p  < 
1.  Thus,  the  effectiveness  of  radius  modulation  is  not  affected  by  momentum 

spread. 

RF  Coupling  Between  Cavities 


We  have  for  simplicity  assumed  that  the  cavities  are  isolated  from  one 
another  electromagnetically  when  the  beam  is  absent.  There  are  several  ways  to 
reduce  the  coupling  among  neighboring  cavities.  The  inductive  coupling  at  the 
annular  slots,  through  which  the  driver  beam  passes,  may  be  cancelled  by  the 
capacitive  coupling  at  the  center  hole,  and  if  necessary,  by  introducing  additiona 
holes  near  the  rf  electric  field  maximum  (so  as  to  increase  the  capacitive  coupling) 
that  is  close  to  the  outer  wall  of  the  cavity  [Fig.  1].  Alternatively,  conducting 
wires  may  be  inserted  radially  across  the  annular  gap  to  reduce  the  inductive 
coupling.  Multiple  pencil  beams  may  also  be  used  as  the  driver.  These  penci 
beams  pass  through  holes  that  are  distributed  annularly.  In  the  event  that  the 
neighboring  cavities  are  not  completely  isolated  electromapetically,  a  traveling 
wave  formulation  would  be  required;  but  the  radius  modulation  that  is  proposed  in 
this  paper  still  provides  an  external  control  to  ensure  phase  focusing. 

The  presence  of  intense  space  charge  in  the  driver  beam  complicates  matters 
substantially,  as  it  is  known  to  modify  the  rf  characteristics  in  an  unpredictable 
manner.  Such  modifications  include  a  detune  of  the  structure  frequency  and 
modification  of  the  gap  transit-time  factor,  especially  if  a  virtual  cathode  is  on  the 
verge  of  being  formed  [2,  11].  Other  modes  may  be  excited.  Indeed,  mode 
competition  is  a  major  area  that  requires  considerable  attention  in  the  twobetron 

concept. 

Effect  of  Finite  Thickness  in  the  Driver  Beam 


The  effect  of  finite  thickness,  x,  in  the  driver  beam  is  found  to  be  much  less 
important  than  its  mean  radius,  ro.  The  finite  beam  thickness  modifies  Eq.  3  to 
read 


5  --1.249(r^ 


-a)/a-0.05(x/a)^. 


(15) 
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The  last  term  in  Eq.  15  is  usually  much  smaller  than  the  first  term  even  if  x  and 
the  radius  modulation  amplitude.  A,  are  of  the  same  order  of  magnitude  (Eq.  6). 

CONCLUDING  REMARKS 

In  summary,  the  twobetron  has  the  potential  of  converting  many  existing 
pulse  power  systems  into  compact  rf  accelerators  that  are  suitable  for  industrial 
and  medical  applications.  The  driver  beam  is  a  modulated  electron  beam  of 
annular  shape  and  low  energy.  The  secondary  beam  is  an  on-axis  pencil  beam. 
The  secondary  beam  may  reach  an  energy  up  to  10  MeV  in  one  to  two  meters. 
Transformer  ratio  on  the  order  of  ten  is  considered  feasible  for  each  stage.  Phase 
focusing  and  energy  tunability  of  the  accelerated  beam  may  be  provided  by  an 
external  magnetic  field,  which  controls  the  radius  of  the  primary  beam. 

Excitation  of  the  undesirable  modes  by  the  driver  beam  is  perhaps  the  single 
most  important  issue  in  the  twobetron  concept. 
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It  is  shown  that  an  annular  electron  beam  may  carry  six  times  as  much  current  as  a  pencil  beam 
for  the  same  beam  breakup  (BBU)  growth.  This  finding  suggests  that  the  rf  magnetic  field  of  the 
breakup  mode  is  far  more  important  than  the  rf  electric  field  in  the  excitation  of  BBU. 
proof-of-principle  experiment  is  suggested,  and  the  implications  explored. 


Annular  electron  beams  have  the  capability  of  carrying 
a  much  higher  current  than  a  pencil  beam.  Besides  the 
obvious  fact  that  annular  beams  have  a  larger  cross- 
sectional  area,  their  limiting  currents  are  significantly 
higher  than  those  of  a  pencil  beam  when  placed  in  a  me¬ 
tallic  drift  tube.  For  this  and  other  reasons,  annular  beams 
have  recently  been  chosen  as  the  preferred  geometry  to 
generate  coherent,  ultrahigh  power  microwaves.  ’  They 
have  also  been  used  as  the  primary  beam  in  several  two- 
beam  accelerator”  configurations.^’"^  These  annular  beams 
either  encounter  a  sequence  of  modulating  gaps,  or  simply 
glaze  by  a  slow  wave  structure  to  generate  a  wake  field  in 
the  case  of  two-beam  accelerators.^  The  beam  radius,  the 
pill  box  radius,  and  the  slow  wave  structure  radius  may  all 
be  of  the  same  order  of  magnitude.  The  high  current  would 
then  lead  to  the  beam  breakup  instability  (BBU)^  and 
this  concern  motivates  the  present  study. 

BBU  is  usually  analyzed  for  a  pencil  beam  propagating 
along  the  center  axis  of  a  sequence  of  accelerating  cavities. 
Many  BBU  calculations  of  practical  interest  assume  that 
the  accelerating  unit  is  the  familiar  cylindrical  pillbox  cav¬ 
ity  and  that  the  dominant  deflecting  mode  is  the  TM^o 
mode.^’^’^^  Extension  to  an  annular  beam  is  straightfor¬ 
ward.  Nevertheless,  this  calculation  leads  to  several  unex¬ 
pected  results  and  provides  some  new  insights  into  BBU,  to 
be  reported  in  this  communication. 

It  is  well  known  that  BBU  is  excited  by  the  combined 
action  of  the  rf  magnetic  field  (Bj)  and  the  rf  electric  field 
(Ej)  of  the  deflecting  modes:^  causes  beam  deflection 
through  the  Lorentz  force  and  Ej  causes  mode  amplifica¬ 
tion  through  the  work  done  on  the  mode  by  the  beam 
current  J.  Our  calculation  strongly  suggests  that  Bj  is 
much  more  critical  than  Ei  in  contributing  to  BBU 
growth.  Thus,  an  annular  beam  strategically  placed  near 
the  minimum  of  the  rf  magnetic  field  would  suffer  far  less 
beam  breakup  growth  than  a  pencil  beam  that  is  centered 
on  the  cavity  axis,  where  the  magnetic  field  is  large  and  the 
axial  electric  field  is  small.  By  the  same  argument,  placing 
the  annular  beam  very  close  to  the  wall  of  a  metallic  drift 
tube,  at  which  the  axial  electric  field  is  vanishingly  small, 
cannot  eliminate  BBU  growth  because  of  the  substantial 
deflecting  magnetic  field  generated  by  the  wall  current. 
Toward  the  end  of  this  communication,  we  propose  an 
experiment  which  would  unambiguously  test  the  relative 
importance  between  the  rf  magnetic  field  and  the  rf  axial 
electric  field,  as  discussed  here. 


Consider  an  infinitesimally  thin  annular  beam  of  radius 
/•q  inside  a  cylindrical  pillbox  of  radius  6.  The  beam  carries 
a  total  current  /  and  coasts  at  velocity  Vq  corre¬ 

sponding  relativistic  factors  y  and  p.  The  drift  tube  is 
loaded  with  a  slow  wave  structure,  modeled  by  a  series  of 
cylindrical  pillbox  cavities,  each  of  which  supports  the 
nonaxisymmetric  TMjio  mode.^'^’^’^^  The  interaction  be¬ 
tween  this  mode  and  the  beam  causes  BBU  to  be  excited. 
In  the  limit  rQ-^O,  this  is  the  basic  model  of  BBU  for  a 
pencil  beam.  Since  we  are  comparing  the  strength  of  BBU 
interaction  for  different  values  of  ro,  we  pretend  that  mag¬ 
netic  focusing  is  absent  and  that  the  quality  factor  Q  of  the 
deflecting  mode  is  infinite. 

Let  Ai=  zq{t)  {cos  d)E(r)  be  the  vector  potential  of 
the  deflecting  dipole  mode  in  a  cavity.  For  the  fundamental 
TMno  niode,  E{r)=Ji{pr)  represents  the  radial  depen¬ 
dence  of  the  axial  electric  field  with  Jy  being  the  Bessel 
function  of  order  one  and  p=3,S32/b,  The  corresponding 
magnetic  field  is  Bj  — VXAj.  The  action  of  this  mode  on 
the  beam  is  calculated  as  follows. 

We  divide  the  annular  beam  into  N  azimuthal  seg¬ 
ments  (A^  large).  The  ith  segment  is  located  at  r=ro, 
O—Bi—l'tri/N,  in  the  unperturbed  state  but  is  displaced 
radially  by  and  azimuthally  by  when  the  deflecting 
mode  is  present.  The  linearized  force  law  yields 

-y{co-kvQ)^^i={e/mo){vo/c)qE'{ro)cos  O,-,  (1) 

e/m^ ) 

X  {vQ/c)q[E{rQ)/rQ\s\n  0/,  (2) 

where  the  right-hand  sides  represent  the  components  of  the 
Lorentz  force  that  causes  beam  deflection.  In  writing  Eqs. 
(1)  and  (2),  we  have  assumed  a  wave-like  solution 
cxy\j{(ot—kz)]  for  the  disturbances,  with  —  1,  and  we 
have  used  a  prime  to  denote  derivative  with  respect  to  the 
argument. 

The  instantaneous  current  J  on  the  iih  current  fila¬ 
ment  is 

J,irj)=z--^-5ir-ro-^i)8(^9-ei-y^,  (3) 

where  8  is  the  Dirac  delta  function.  The  work  done  by  this 
current  filament  on  the  deflecting  mode  is  proportional  to 

1PF,=  J  dVA^-Ji,  (4) 
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FIG-  1.  Comparison  of  the  BBU  coupling  constant  e  between  an  annular 
beam  of  radius  tq  and  an  on-axis  p>encil  beam  (ro-*0)  with  the  same  total 
current. 


where  the  volume  integral  is  performed  over  the  cavity.  In 
evaluating  IV (i  we  should  retain  only  the  rf  component  of  J/ 
in  Eq.  (3),  since  only  the  rf  current  performs  work  on  the 
breakup  mode.  Upon  substituting  Eqs.  (1)— (3)  into  Eq. 
(4),  and  summing  over  all  /,  we  find  the  total  work  done 


N 


w=  X  Wi= 


Lie  Vq  [£''(/b) 
^0  ^  ^ 


(5) 


apart  from  a  multiplicative  constant  that  is  independent  of 
the  beam’s  equilibrium  position  Kq.  This  energy  transfer 
leads  to  growth  of  the  BBU  mode,  which  is  described  by 
the  BBU  dispersion  relation:^ 


kvo)^  =  —  IcoQe  —  —  Ico^e^^e/  Fq)  ,  (6) 


where  e  is  the  coupling  constant  and  coq  is  the  breakup 
mode  frequency.  In  writing  the  last  form  of  Eq.  (6),  we 
normalize  e  in  terms  of  6o>  coupling  constant  for  an 
on-axis,  pencil  beam  (rQ— ►O).  For  the  TMhq  mode, 
E^Jy{pr)  and  6o=OA22ip/r)(I/l  kA).  It  is  clear  from 
Eq.  (5)  that 


e  1 

which  compares  the  BBU  strength  between  an  annular 
beam  and  a  pencil  beam  of  the  same  current.  Note  that  this 
ratio  reduces  to  unity  in  the  limit 

Equation  ( 7 )  is  plotted  in  Fig.  1  as  a  function  of  f'c/b. 
It  is  seen  from  this  figure  that  €/€q  may  be  as  small  as  0.17 
when  the  annular  beam  is  located  at  ro=0.56^.  Note  also 
that  this  location  coincides  with  the  minimum  of  the  rf 
magnetic  field  of  the  deflecting  mode.  What  this  means  is 
that  an  annular  beam  placed  at  this  location  can  carry  as 
much  as  1/0.17  =  6  times  the  current  as  an  on-axis  pencil 
beam,  and  suffer  the  same  BBU  growth.  Another  point 
worth  noting  is  that  BBU  growth  retains  significant 
strength  even  if  the  annular  beam  is  very  close  to  the  wall 
of  the  drift  tube  (cf.  rQ-^bin  Fig.  1).  This  result  is  unex¬ 


pected  since  Ej— >0  near  a  metallic  wall.  As  a  result, 
and,  superficially,  one  could  hardly  expect  any 
transfer  of  power  from  the  beam  to  drive  the  breakup 
mode.^  The  finite  BBU  strength  ro^b  is  another  strong 
indication  that  the  deflecting  magnetic  field  is  far  more 
important  than  the  axial  rf  electric  field  in  driving  BBU. 

The  importance  of  the  rf  magnetic  field  can  be  tested  in 
an  experiment  in  which  a  pencil  beam  is  focused  by  a 
solenoidal  magnetic  field  and  is  made  to  pass  through  a 
sequence  of  pillbox  cavities,  in  which  the  first  cavity  is 
primed  with  microwaves  at  the  TM^o  mode.  BBU 
growth  is  monitored  at  the  last  cavity,  before  the  beam 
exit.  The  above  theory  then  predicts  the  unusual  feature 
that  BBU  growth  should  be  much  less  if  the  pencil  beam  is 
placed  off-axis,  than  if  the  pencil  beam  were  on-axis.  The 
BBU  growth  should  be  minimum  if  this  pencil  beam  is 
placed  at  a  distance  of  about  0.56  of  the  pillbox  radius, 
where  the  rf  magnetic  field  is  minimum. 

We  also  repeated  the  calculations  for  the  higher  order 
radial  modes:  TMi2o»  TM130,  TM140,  and  TMjso*  Fixing 
ro/6  =  0.56,  the  ratio  e/e^  equals  0.16,  0.012,  0.037,  and 
0.013  for  these  four  higher  order  modes,  respectively. 
Thus,  the  annular  beam  still  suffers  substantially  lower 
BBU  growth,  in  the  higher  order  deflecting  modes,  than  an 
on-axis  pencil  beam  of  the  same  current. 

In  conclusion,  the  rf  magnetic  field  is  found  to  be  much 
more  important  than  the  rf  electric  field  in  contributing  to 
BBU  growth.  A  simple  proof-of-principle  experiment  is 
proposed  to  test  this  new  finding.  Annular  beams  are  far 
more  stable  than  an  on-axis  pencil  beam,  as  a  result. 
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^The  rate  of  power  transfer  is  proportional  to  the  gradient  of  Ej,  rather 
than  El  itself.  This  may  be  seen  when  one  substitutes  Eq.  (3)  into  Eq. 
(4)  and  performs  integration  by  parts. 
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"BBU  growth  on  a  pencil  beam  that  is  placed  off-center  can  be  easily 
calculated  by  using  Eq.  (4)  instead  of  Eq.  (5).  We  pretend  that  the 
total  beam  current  is  carried  by  the  /th  filament  that  enters  Eq.  (4). 
Although  the  BBU  growth  of  such  an  off-center  beam  depends  on  6^  its 
coupling  constant  €  is  still  much  less  than  €;q»  value  for  an  on-axis 
beam. 
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Abstract 


By  the  use  of  a  simple  model,  we  explicitly  incorporate  the  coupling  between  the  driver 
cavity  and  the  booster  cavity  in  a  relativistic  klystron  amplifier  (RKA).  We  show  that  this  RKA 
configuration  may  turn  into  an  injection  locked  oscillator  only  when  the  beam  current  is  sufficiently 
high.  Other  features  revealed  by  this  model  include:  the  downshifted  frequency  mode  ("0"  mode) 
is  unstable  whereas  the  upshifted  frequency  mode  ("tc"  mode)  is  stable;  the  growth  rate  of  the  0 
mode  is  relatively  mild  so  that  the  oscillation  can  start  only  in  an  injection  locked  mode;  the 
oscillation  does  not  require  the  presence  of  reflected  electrons;  and  the  separation  of  the  cavities 
must  be  sufficiently  short.  These,  and  other  features,  are  found  to  be  in  qualitative  agreement  with 
the  recent  experiments  on  the  injection  locked  relativistic  klystron  oscillator  (RKO)  that  were 
conducted  at  the  Phillips  Laboratory. 
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T.  Introduction 


Current  modulation  of  an  intense  relativistic  electron  beam  (lEEB)  remains  a  challenging 
and  active  research  problem  [1]  with  applications  to  high  power  microwave  (HPM)  production  [2- 
4],  and  to  some  accelerator  schemes  that  use  an  IREB  as  a  driver  [5,6].  Recently,  a  series  of 
experiments  were  performed  at  Phillips  Laboratory  [7],  where  an  IREB  was  placed  in  a  relaUvistic 
klystron  amplifier  (RKA)  geometry  (2.4),  with  the  driver  and  booster  caviues  sufficienUy  close  to 
couple  to  one  another.  This  cavity  coupling  caused  the  smicture  to  oscillate  at  a  frequency  different 
from  the  driver  frequency.  This  oscillaUon  persisted  after  the  external  drive  was  shut  off,  but  its 
appearance  has  always  required  an  input  drive  and  a  sufficienUy  large  DC  current.  Thus,  the 
device  operated  as  an  injection-locked  relativistic  klystron  oscillator  (RKO).  It  provided  40%  to 
60%  of  current  modulation  on  the  DC  beam  when  the  gap  separation  is  on  the  order  of  only  10  cm, 
as  opposed  to  30  cm  in  a  typical  RKA  geometry  with  similar  frequency  (at  1.3  GHz),  beam  current 
and  drift  tube  diameter.  Another  interesting  feature  is  that  virtual  cathodes  did  not  seem  to  have 
been  foimed  for  the  operation  of  this  RKO.  The  feedback,  therefore,  was  due  only  to  the 
electromagnetic  coupling  among  the  two  cavittes,  and  not  to  the  reflected  electrons  as  in  previous 

works  [8]. 

In  this  paper,  we  present  a  simple  analytic  model,  explicitly  including  the  effect  of  cavity 
coupling.  The  results  of  this  model  demonstrate  several  features  that  were  observed  m  the  RKO 
experiments  mentioned  above.  They  are  summarized  in  the  abstract  and  are  discussed  further 


below. 


TT.  Model 


Consider  an  annular  intense  relativistic  electron  beam  (IREB)  with  radius  rb  propagating  in 
a  drift  tube  of  radius  rw  and  passing  by  two  cavities  [Fig.  1].  The  first  cavity,  with  gap  voltage 
A(t),  is  driven  by  a^i  external  source.  The  second  cavity,  with  gap  voltage  B(t)  and  located  at  a 
distance  d  downstream,  is  driven  by  the  beam.  When  d  is  sufficiently  small,  there  would  be 
coupling  between  A  and  B  even  if  the  drift  tube  is  cutoff  to  the  frequency  of  operation.  In  the 
absence  of  the  beam,  the  steady  state  voltages  A  and  B  are  either  in  phase  ("0"  mode)  or  180°  out 
of  phase  ("pi"  mode).  That  is,  there  is  no  phase  delay  between  the  two  cavities,  when  the  beam  is 
absent,  because  the  drift  tube  is  cutoff  [9]. 

To  study  this  effect,  we  start  by  writing  down  the  circuit  equations  of  the  two  cavities,  each 
having  a  natural  frequency  C0o>  and  quality  factor  Q.  The  evolution  of  the  gap  voltages  A  and  B  is 

governed  by 


cOq  d 

dt^  Q  dt 


■to; 


A(t)  =  O)oCB(t)  , 


(1) 


'd^ 

1-^A+ 

dt^' 

Q  dt 

B(t)  =  +  OL)oCA(t) 


(2) 


where  C  («  1)  is  the  dimensionless,  real  constant  measuring  the  degree  of  coupling  between  the 
two  gaps  when  the  beam  is  absent  [9],  and  Z  is  defined  such  that  ZQ  is  the  impedance  (in  Ohms)  at 
the  second  gap.  Equation  (1)  expresses  the  excitation  of  the  gap  voltage  A  by  the  gap  voltage  B  as 
a  result  of  cavity  coupling  C.  Since  we  are  concentrating  on  the  RKO  operation,  we  ignore  the 
external  rf  drive  on  the  gap  voltage  A,  and  we  envision  the  effect  of  this  external  drive  to  enter  only 
as  the  initial  condition  on  A,  at  the  instant  when  the  external  drive  is  shut  off.  Equation  (2) 


describes  the  excitation  of  the  second  gap  by  the  rf  current,  I,  and  by  the  first  gap  voltage  A 
because  of  the  cavity  coupling. 


As  in  the  usual  klystron  theory,  the  rf  current  I  at  the  second  gap  is  due  to  the  voltage  at  the 
first  gap.  In  the  case  of  an  IREB  of  DC  current  Iq,  it  is  given  by  [1] 

In  Eq.  (3),  R  is  the  rf  beam  impedance  (in  Ohms),  0  =  cod/pc  is  the  phase  delay  of  the  beam  mode, 
and 
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R  = 


SllkV 
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a  = 


Ip 

IsPy' 


I3  =8.5kA/£n(rw /Pb)  > 


(4) 


where  p  and  y  are  the  usual  relativistic  mass  factors  associated  with  the  DC  beam,  and  c  is  the 
speed  of  light.  Upon  inserting  Eq,  (3)  into  Eq.  (2),  we  obtain  two  equations  (1).  (2)  in  two 
unknowns  A,  B.  Assuming  exp(j®t)  dependence  for  both  A  and  B,  we  obtain  the  dispersion 

relation  for  o)  which  may  easily  be  solved  to  yield 


(0  =  ® 


1  +  J-  ±  £1  +  — sin(kpd)e 
2Q  2  V  CR 


-je 


(5) 


Equation  (5)  gives  the  temporal  growth  rate  of  the  coupled-cavity  RKO  in  terms  of  the  circuit 
parameters  (cOp,  Q,  C,  d,  Z)  and  the  beam  parameters  (R,  kp,  0).  In  the  next  section,  we  present 
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the  numerical  results  using  parameters  similar  to  those  in  the  experiments  [7].  We  shall  also 
describe  the  interesting  features  revealed  by  the  dispersion  relation  (5). 


TTI.  Numerical  Results 

We  shall  first  establish  the  numerical  values  of  the  parameters  from  the  experiments  [7J.  In 
the  cold  tests,  it  was  observed  that  when  the  two  cavities  are  coupled,  there  was  a  frequency  shift 
of  ±7  MHz  from  the  natural  frequency  of  1.27  GHz  that  was  measured  for  an  isolated  cavity. 
Thus,  when  the  beam  is  absent,  R  -o  ~  from  Eq.  (4)  and  we  obtain  from  Eq.  (5)  the  cold  tube 

coupled-cavity  natural  frequencies 


with  cOo  =  27C  X  1.27  GHz  and  cOoC/2  =  27t  x  7  MHz.  This  yields  the  coupling  constant  C  = 

0.01 12.  In  Eq.  (6),  the  (+)  sign  corresponds  to  the  mode  and  the  (-)  sign  corresponds  to  the 
"0"  mode.  Substitution  of  Eq.  (6)  into  Eq.  (1),  in  which  d/dt  is  replaced  by  j(0,  shows  that  the  gap 
voltages  A  and  B  are  180®  out  of  phase  for  the  "n"  mode,  and  that  they  are  in  phase  for  the  "0" 

mode. 

The  value  of  Q  to  be  used  is  uncertain.  There  is  the  cold  tube  value,  in  excess  of  400 
according  to  measurements  [7],  and  there  is  the  substantially  reduced  value  of  the  beam-loaded  Q, 
which  depends  on  the  beam  current.  Let  us  arbitrarily  set  the  beam-loaded  Q  to  be  100,  say.  The 
value  of  Z  may  be  estimated  from  Eq.  (2)  which  shows  that,  in  the  absence  of  cavity  coupling 
[i.e.,  C  =  0],  the  magnitude  of  the  gap  voltage  B  equals  to  QZ  times  the  magnitude  of  the  rf  current 
I  there.  From  this  definition  of  Z,  we  estimate  that  Z  =  1.6  Q  when  Q  is  assigned  to  be  100.  The 
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beam  parameters  R,  kp,  and  6,  ate  deteimined  as  follows.  We  set  the  beam  voltage  at  400  keV, 
beam  current  I„  =  12  kA,  with  mean  beam  radius  rp  =  6.85  cm  in  a  drift  tube  of  radius  r*  =  7.65 
cm.  Using  these  parameters,  we  find  Y  =  1 .783,  P  =  0.828, 1,  =  77  kA.  a  =  0.033,  kp  =  0.0398 
cm-l  and  R  =  20.4  £2.  If  we  set  the  gap  separation  at  d  =  11  cm,  then  kpd  =  0.437  and  0  =  3.55. 

Using  these  parameters  in  Eq.  (5),  we  find 


®  =  Q)o(l-0023  +  0.0127j) 
a)  =  ©o(0-9976-0.00266j) 


"7t"  mode 
"0"  mode 


(7a) 

(7b) 


Equations  (7a)  and  (7b)  show  the  following. 


A.  The  "k"  mode  is  stable  and  the  "0"  mode  is  unstable. 

B.  For  the  "0"  mode,  there  is  a  downshift  in  frequency  from  fo  =  (oJ2k  =  1.27  GHz,  by 

the  amount  of  1.27  GHz  x  (1  -  0.9976)  =  3.05  MHz.  That  is,  the  mode  is  upshifted  from  the 
cold-tube  "0"  mode  frequency  by  the  amount  of  7  MHz  -  3.05  MHz  =  3.  95  MHz,  and  this 
upshift  is  beam  induced. 

C.  The  growth  of  the  "0"  mode,  according  to  Eq.  (7b),  is  relatively  mild.  The  total  number  of 

e-folds  in  a  time  x  =  100  ns  is  (DiX  =  2.12.  This  mild  growth  implies  that,  for  an  IREB  whose 
pulselength  is  on  the  order  of  100  ns,  the  manifestation  of  RKO  behavior  may  require  an 
external  drive  so  that  there  is  already  a  significant  rf  gap  voltage  by  the  time  this  external  rf 

drive  is  shut  off. 

D.  The  "0"  mode  is  unstable  only  if  the  DC  beam  current  is  sufficiently  high.  This  is  obvious 

since  in  the  limit  of  zero  current,  the  mode  is  damped  [cf.  Eq.  (6)].  Using  the  parameters 
given  in  the  paragraph  preceeding  Eq.  (7a),  we  find  that  the  threshold  DC  beam  current  for  the 
onset  of  growth  for  the  "0"  mode  is  7.2  kA.  Similar  levels  of  threshold  current  were  observed 


in  the  experiments. 


E.  Upon  inserting  Eq.  (7b)  into  Eq.  (1),  in  which  d/dt  is  replaced  by  jffi,  we  find  the  ratio  of  the 
gap  voltage  to  be  |B  /  A]  =  1.43  using  the  parameters  in  the  paragraph  preceedmg  Eq  (7a). 

This  number  is  in  reasonable  agreement  with  experimental  observations. 

Points  A  -  E  are  features  revealed  by  the  experiments  [7].  We  wish  to  add  the  following 

points  which  also  seem  to  corroborate  the  experimental  results. 

F.  The  instability  does  not  require  the  formation  of  a  virtual  cathode.  That  is,  there  is  no  need  to 
invoke  reflected  electrons  to  provide  the  feedback  that  is  usually  required  for  oscillation  [8]. 

G.  The  cavity  separation  must  be  sufficiently  short  to  provide  appreciable  couphng  between  the 

two  cavities. 

H.  The  operation  of  this  injection-locked  RKO  is  restricted  only  to  a  narrow  range  of  gap 
separation,  d,  once  the  other  circuit  parameters,  beam  parameters,  and  the  level  of  rf  drive  are 
fixed.  There  are  three  reasons  for  this:  (i)  The  cavity  coupling  diminishes  rapidly  with 
increasing  d,  as  the  drift  tube  is  below  cutoff,  (ii)  If  d  is  too  small,  there  is  little  buildup  in  the 
current  modulation  from  the  driver  cavity  to  the  second  one  downstream,  (iii)  The  ballistic 
phase  e  that  appears  in  Eq.  (5)  is  proportional  to  d.  Note  that  (i)  is  quantified  in  Eq.  (5) 
through  the  coupling  constant  C  and  (ii)  is  quantified  in  Eq.  (5)  through  the  factor  sin(kpd). 


TV.  Concluding  Remarks. 


While  we  are  able  to  construct  a  relatively  simple  model  which  appears  to  be  capable  of 
explaining  many  features  observed  in  the  recent  injection-locked  RKO  experiments,  it  must  be  kept 
in  mind  that  there  are  many  uncertainties  in  the  model.  Chief  among  them  are  the  values  of  Q  and 
the  coupling  constant  C.  The  beam-loaded  Q  depends  on  the  current  modulation,  which  depends 
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not  only  on  the  DC  current,  but  also  on  the  level  of  rf  drive  and  the  gap  separation  d.  The  value  of 
C  is  a  sensitive  function  of  gap  separation,  and  the  one  we  use  in  this  paper  is  inferred  from  the 
cold  tube  measurements.  There  is  also  uncertainty  in  the  value  of  the  gap  impedence  Z  because  the 
"gap  transit  time  factor",  which  is  significantly  affected  by  the  DC  space  charge  effects  associated 
with  an  IREB,  will  make  the  determination  of  Z  far  from  a  trivial  matter  [10].  Unfortunately,  the 
conceptually  simnle  mechanism  proposed  in  this  paper  may  not  be  easy  to  verify  in  particle-in-cell 
codes.  The  main  reason  is  that  the  mild  growth  envisioned  would  require  a  long  simulation  time. 
To  shorten  the  simulation  time  by  raising  the  beam  current  runs  into  the  possibility  of  triggenng 
virtual  cathode  formation.  In  addition,  there  may  be  marked  difference  between  the  "numerical  Q", 
the  cold-tube  Q,  and  the  hot-tube  Q  in  the  experiments.  As  explained  above,  our  model  shows  that 
the  threshold  current  may  depend  sensitively  both  on  Q  and  on  the  gap  separation  d. 

In  spite  of  the  great  uncertainties  in  several  crucial  parameters,  we  find  it  remarkable  that  a 
reasonable  choice  of  parameters  does  yield  reasonable  agreement  with  observations,  based  just  on 
the  simple  analytic  model. 
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